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The Standard Model (SM) of particle physics contains about two dozen parameters ~ such as particle 
masses whose origins are still unknown and cannot be predicted, but whose values are constrained 
through their interactions. In particular, the masses of the top (t) quark (Aft) and W boson (Mw) 
[V\ constrain the mass of the long-hypothesized, but thus far not observed, Higgs boson. A precise 
measurement of the top-quark mass can therefore point to where to look for the Higgs, and indeed 
whether the hypothesis of a SM Higgs is consistent with experimental data. Since top quarks are 
produced in pairs and decay in only ~ 10~^^ s into various final states, reconstructing their mass 
from their decay products is very challenging. Here we report a technique that extracts far more 
(^ information from each top-quark event and yields a greatly improved precision on the top mass of 

f^ ±5.3 GeV/c^, compared to previous measurements [2j. When our new result is combined with our 

PsJ published measurement in a complementary decay mode j^ and with the only other measurements 

^ available [2j, the new world average for Mt becomes 178.0 ±4.3 GeV/c^ \4\. As a result, the most likely 

^ Higgs mass increases from the experimentally excluded [5j value of 96 GeV/c^ [6\ to 117 GeV/c^, which 

^~~i is beyond current experimental sensitivity. The upper limit on the Higgs mass at 95% confidence level 
,_! ; is raised from 219 GeV/c^ to 251 GeV/c^. 

I The discovery of the top quark in 1995 served as one of the major confirmations of the vahdity of the SM |3jl3l- Of 

I . its many parameters, the mass of the top quark, in particular, reflects some of the most crucial aspects of the Model. 

^ ■ This is because, in principle, the top quark is point-hkc and should be massless; yet, through its interactions with the 
,— I , hypothesized Higgs field, the physical mass of the top quark appears to be about the mass of a gold nucleus. Because 
C^ ' it is so heavy, the top quark (along with the W boson) provides an unusually sensitive tool for investigating the Higgs 
^D ' field. Mw is known to a precision of 0.05%, while the uncertainty on Mt is at the 3% level 1]. Improvements in 
both measurements are required to restrict further the allowed range of mass for the Higgs; however, given the large 
uncertainty in the top mass, an improvement in its precision is particularly important. As has been pointed out 
recently, e.g., in Refs. |3,[i3) a potential problem for the SM is that, based on the presently accepted mass of the top 
quark, the most likely value of the Higgs mass 6i] lies in a range that has already been excluded by experiment |a| . 

^ . Precise knowledge of the Higgs mass is crucial for our understanding of the SM and any possible new physics beyond 
it. For example, in a large class of supersymmetric models (theoretically preferred solutions to the deficiencies of the 

Qh! SM), the Higgs mass has to be less than w 135 GeV/c^. While, unlike the SM, supersymmetry predicts more than 

^5 one Higgs boson, the properties of the lightest one are expected to be essentially the same as those for the SM Higgs 

boson. Thus, if the SM-like Higgs is heavier than « 135 GeV/c^, it would disfavor a large class of supersymmetric 
models. In addition, the current limits on supersymmetric particles from LEP [llj are extremely sensitive to the 
mass of the top quark. In fact, for a top-quark mass greater than 179 GeV/c^, the bounds on one of the major 
supersymmetry parameters, tan/3, which relates the properties of the SM-like Higgs boson and its heavier partners, 

C^ ' would disappear completely |l2l |. Hence, in addition to the impact on searches for the Higgs boson, other important 
consequences call for improved precision on the top mass, and this goal is the main subject of this paper. 

The D0 experiment at the Fermilab Tevatron has studied a sample of ti events produced in proton-antiproton 
(pp) interactions [13. The total energy of 1.8 TeV released in a head-on collision of a 900 GeV p and 900 GeV p 
is almost as large as the rest energy of ten gold nuclei. Each top (antitop) quark decays almost immediately into 
a bottom b (b) quark and a W^ {W~) boson, and we have reexamined those events in which one of the W bosons 
decays into a charged lepton (electron or muon) and a neutrino, and the other W into a quark and an antiquark 
(see Figure . These events and their selection criteria are identical to those used to extract the mass of the top 
quark in our previous publication, and correspond to an integrated luminosity of 125 events/pb. (That is, given 
the production cross section of the tt in pp collisions at 1.8 TeV of 5.7 pb, as measured by D0 [lj|j these data 
correspond to approximately 700 produced tt pairs, a fraction of which is fully detected in various possible decay 
modes. Approximately 30% of these correspond to the "lepton + jets" topology categorized in Figure [21 where "jet" 
refers to products of the fragmentation of a quark into a coUimated group of particles that are emitted along the 
quark's original direction.) The main background processes correspond to multijet production (20%), where one of 
the jets is reconstructed incorrectly as a lepton, and the W+ jets production with Icptonic W decays (80%), which 
has the same topology as the tt signal. 

The previous D0 top-quark mass measurement in this lepton + jets channel is Mt = 173.3 ± 5.6 (stat) ± 
5.5 (syst) GeV/c^, and is based on 91 candidate events. Information pertaining to the detector and to the older 
analysis can be found in Refs. [13 and [13, respectively. 
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The new method of the top-quark mass measurement is similar to one suggested previously [la, uM fo'" ^^ dilepton 
decay channels (where both W bosons decay leptonically), and used in previous mass analyses of dilepton events 0], 
and akin to an approach suggested for the measurement of the mass of the W boson at LEP [ll, llS l20| ■ The critical 
differences from previous analyses in the lepton + jets decay channel lie in: (i) the assignment of more weight to 
events that are well measured or more likely to correspond to tt signal, and (ii) the handling of the combinations of 
final-state objects (lepton, jets, and imbalance in transverse momentum, the latter being a signature for an undetected 
neutrino) and their identification with top-quark decay products in an event (e.g., from ambiguity in choosing jets 
that correspond to 6 or 6 quarks from the decays of the t and i quarks). Also, since leading-order matrix elements 
were used to calculate the event weights, only events with exactly four jets are kept in this analysis, resulting in a 
candidate sample of 71 events. Although we are left with fewer events, the new method for extracting the mass of 
the top quark provides substantial improvement in both statistical and systematic uncertainties. 

We calculate as a function of top mass the differential probability that the measured variables in any event corre- 
spond to signal. The maximum in the product of these individual event probabilities provides the best estimate of the 
mass of the top quark in the data sample. The impact of biases from imperfections in the detector and event recon- 
struction algorithms is taken into account in two ways. Geometric acceptance, trigger efficiencies, event selection, etc., 
enter through a multiplicative acceptance function that is independent of Mt . Because the angular directions of all 
the objects in the event, as well as the electron momentum, are measured with high precision, their measured values 
are used directly in the calculation of the probability that any event corresponds to tt or background production. The 
known momentum resolution is used to account for uncertainties in measurements of jet energies and muon momenta. 

As in the previous analysis [13, momentum conservation in 7+jet events is used to check that the energies of jets 
in the experiment agree with Monte Carlo (MC) simulation. This calibration has an uncertainty SE = (0.025E' -I- 0.5 
GeV). Consequently, all jet energies in our sample are rescaled by ±5E, the analysis redone, and half of the difference 
in the two rescaled results for Mt (SMt = 3.3 GeV/c^) is taken as a systematic uncertainty from this source. All 
other contributions to systematic uncertainty: Monte Carlo modeling of signal {SMt = 1.1 GeV/c^) and background 
{6Mt = 1.0 GeV/c^), effect of calorimeter noise and event pile-up {SMt = 1.3 GeV/c^), and other corrections from 
top-quark mass extraction {SMt = 0.6 GeV/c^) are much smaller, and discussed in detail in Refs. '23,|22|. It should be 
noted that the new mass measurement method provides a significant (« 40%, from ±5.5 to ±3.9 GeV/c^) reduction 
in systematic uncertainty, which is ultimately dominated by the measurement of jet energies. For details on the new 
analysis, see the Methods section. 

The final result is Alt = 180.1 ± 3.6 (stat) ± 3.9 (syst) GeV/c^. The improvement in statistical uncertainty over 
our previous measurement is equivalent to collecting a factor of 2.4 as much data. Combining the statistical and 
systematic uncertainties in quadrature, we obtain Mt = 180.1 ± 5.3 GeV/c^, which is consistent with our previous 
measurement in the same channel (at w 1.4 standard deviations), and has a precision comparable to all previous 
top-quark mass measurements combined ™. 

The new measurement can be combined with that obtained for the dilepton sample that was also collected at D0 
during Run I (M* = 168.4 ± 12.3 (stat) ± 3.6 (syst) GeV/c^), to yield the new D0 average for the mass of the top 
quark: 

Mt = 179.0 ± 5.1 GeV/c2 (D0). (1) 

Combining this with measurements from the CDF experiment |3, provides a new "world average" (based on all 
measurements available) for the top-quark mass |j]: 

Mt = 178.0 ± 4.3 GeV/c^ (AU available data), (2) 

dominated by our new measurement. This new world average shifts the best- fit value of the expected Higgs mass from 
96 GeV/c^ to 117 GeV/c^ (see Figure O, which is now outside of the experimentally-excluded region, yet accessible in 
the current run of the Tevatron and at future runs at the Large Hadron Collider (LHC), currently under construction 
at CERN. (The upper limit on the Higgs mass at 95% confidence level changes from 219 GeV/c^ to 251 GeV/c^.) 
Figure 121 shows the effect of using only the new D0 top mass for fits to the Higgs mass, and indicates a best value of 
123 GeV/c^ and the upper limit of 277 GeV/c^ at 95% confidence level. It should be noted that the horizontal scale 
in Figure 131 is logarithmic, and the limits on the Higgs boson mass are therefore asymmetric. 

The new method is already being applied to data being collected by the CDF and D0 experiments at the new 
run of the Fermilab Tevatron and should provide even higher precision on the determination of the top-quark mass, 
equivalent to more than a doubling of the data sample, relative to using the conventional method. An ultimate 
precision of at « 2 GeV/c^ on the top-quark mass is expected to be reached in several years of Tevatron operation. 
Further improvement may eventually come from the LHC. 



Methods 

The probability density as a function of Mt can be written as a convolution of the calculable cross section and any 
effects from measurement resolution: 

Fix, Ah) - — |— fda{y,Mt)dqidq2f{qi)fiq2)W{y,x) (3) 

where W{y, x), our general transfer function, is the normalized probability for the measured set of variables x to arise 
from a set of nascent (partonic) variables y, da{y,Mt) is the partonic theoretical differential cross section, f{q) are 
parton distribution functions that reflect the probability of finding any specific interacting quark (antiquark) with 
momentum q within the proton (antiproton), and a{Mt) is the total cross section for producing ti. The integral in 
Eq. (PJ sums over all possible parton states leading to what is observed in the detector. 

The acceptance of the detector is given in terms of a function A{x) that relates the probability Pjn{x,Mt) of 
measuring the observed variables x to their production probability P{x,Mt): Pm{x,Mt) = A{x)P{x,Mt). Effects 
from energy resolution, etc., are taken into account in the transfer function W{y,x). The integrations in Eq. ^ over 
the eleven well-measured variables (three components of charged-lepton momentum and eight jet angles) and the four 
equations of energy-momentum conservation, leave five integrals that must be performed to obtain the probability 
that any event represents ti (or background) production for some specified value of top mass Mt . 

The probability for a ii interpretation can be written as: 

P^,= J- [d'n E IM*f^^^^i^<i>6W-et.(i?pa..t,i?jct), 

where fl represents a set of five integration variables, Mtt is the leading-order matrix element for ti production 
|24ll2a| . f{qi) and f{q2) are the CTEQ4M parton distribution functions for the incident quarks [23, <&6 is the phase- 
space factor for the 6-object final state, and the sum is over all 12 permutations of the jets and all possible neutrino 
solutions. Wjets(-£'part,^jct) corresponds to a function that maps parton-level energies impart to energies measured in 
the detector -Ejet, and is based on MC studies. A similar expression, using a matrix element for 14^-1- jets production 
(the dominant background source) that is independent of Mt, is used to calculate the probability for a background 
interpretation, Pbkg- 

Studies of samples of HERWIG [23 MC events indicate that the new method is capable of providing almost a 
factor of two reduction in the statistical uncertainty on the extracted Mt- These studies also reveal that there is a 
systematic shift in the extracted Mt that depends on the amount of background there is in the data. To minimize 
this effect, a selection is introduced based on the probability that an event represents background. The specific value 
of the Pbkg cutoff is based on MC studies carried out before applying the method to data, and, for a top mass of 175 
GeV/c^, retains 71% of the signal and 30% of the background. A total of 22 data events out of our 71 candidates 
pass this selection. 

The final likelihood as a function of Mt is written as: 

N 

lnL{Mt) = E HclPuix^,Mt) + C2Pbkg(2;.)] -N A{x) [ciPti{x, Mt) + C2Pbkg(a;)] da;, 

z— 1 

The integration is performed using MC methods. The best value of Mt (when L is at its maximum imax) represents 
the most likely mass of top in the final iV-event sample, and the parameters Ci refiect the amounts of signal and 
background. MC studies show that there is a downward shift of 0.5 GeV/c^ in the extracted mass, and this correction 
is applied to the result. Reasonable changes in the cutoff on Pbkg do not have significant impact on Mt- 

Figure 0] shows the value of L{Mt) / L^^x as a function of Mt for the 22 events that pass all selection criteria, after 
correction for the 0.5 GeV/c^ bias in mass. The likelihood is maximized with respect to the parameters Ci at each 
mass point. The Gaussian fit in the figure yields Mt = 180.1 GeV/c^, with a statistical uncertainty of 5Mt = 3.6 
GeV/c^. The systematic uncertainty, dominated by the measurement of jet energies, as discussed above, amounts 
to 5Mt = 3.9 GeV/c^. When added in quadrature to the statistical uncertainty from the fit, it yields the overall 
uncertainty on the new top-quark mass measurement of ±5.3 GeV/c^. 

We wish to note the great number of contributions made by the late Harry Melanson to the D0 experiment, through 
his steady and inspirational leadership of the physics, reconstruction and algorithm efforts. 
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Figure 1: Feynman diagrams for ft production in pp collisions, with subsequent decays into an electron, neutrino, and quarks. 
Diagram (a) (quark-antiquark production) is dominant, but diagram (b) (gluon fusion) contributes ~ 10% to the cross section. 
This particular final state {ei>udbb) is one of the channels used in the analysis. 
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Figure 2: Relative importance of various tt decay modes. The "lepton + jets" channel used in this analysis corresponds to the 
two offset slices of the pie-chart and amounts to 30% of all the ti decays. 
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Figure 3: Current experimental constraints on the mass of the Higgs boson. The x for a global fit to electroweak data using 
the procedure of Ref. 0, is shown as a function of the Higgs mass. The solid line corresponds to the result for the previous 
world average for the top-quark mass of 174. 3± 5.1 GeV/c^, with the blue band indicating the impact of theoretical uncertainty. 
The dotted line shows the result for the new world-averaged Mt of 178.0 ± 4.3 GeV /c? , while the dashed line corresponds to 
using just the new D0 average of 179.0 ±5.1 GeV/c^. The yellow shaded area on the left indicates the region of Higgs masses 
excluded by experiment (> 114.4 GeV/c^ at the 95% confidence level l^). The improved top mass measurement shifts the 
most likely value of the Higgs mass above the experimentally excluded range. 
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Figure 4: Determination of the top-quark mass using the maximum likelihood method. The points represent the likelihood 
of the fit used to extract the top mass, divided by its maximum value, as a function of the mass of the top quark (after a 
correction for a —0.5 GeV/c mass bias, see text). The solid line shows a Gaussian fit to the points. The maximum likelihood 
corresponds to a mass of 180.1 GeV/c^, which is the new D0 measurement of the top mass in the lepton -I- jets channel. The 
hatched band corresponds to the range of ±1 standard deviation, and indicates the ±3.6 GeV/c^ statistical uncertainty of the 
fit. 



